Biomechanics and Modeling in Mechanobiology

Endothelial cells integrate biomechanical and biochemical cues to regulate vascular function, which can be modeled using cell signaling networks. Yet, the extent to which their signaling behavior is determined by network topology versus reaction-specific parameterization remains unclear. We developed a continuous logic-based ordinary differential equation model of vascular endothelial signaling containing 70 species and 123 regulatory reactions assembled from the literature. Model predictions were evaluated against 78 literature-defined directional responses to increased wall stress, wall shear stress, hypoxia, and TNF. We first tested whether a single global parameter triplet describing reaction weight (w), Hill coefficient (n), and half-maximal activation (EC50) could reproduce these observations. Investigator-selected global parameters reproduced 62 of 78 responses (79.49%). Monte Carlo sampling of 100,000 globally constant parameter sets improved agreement to 69 of 78 responses (88.46%), indicating that network topology alone captured most qualitative signaling behavior but remained sensitive to parameter choice. Allowing reaction-specific parameters increased the best agreement to 75 of 78 responses (96.15%). We then used a genetic algorithm to identify the smallest subset of reactions requiring reaction-specific specification while retaining a common constant-parameter spine for the remainder of the network. Complete agreement (100%) with all 78 literature-defined responses was achieved with reaction-specific parameterization of only 10 of 123 reactions. The optimized model was subsequently used to simulate continuous changes in wall shear stress, revealing pathway-dependent nonlinear responses, including monotonic increases in nitric oxide and pSmad1/5, nonmonotonic responses in endothelin-1 and NF-{kappa}B, and a biphasic pSmad2/3 response. These findings show that most observed endothelial signaling behavior can be reproduced by a fixed network topology with broadly shared reaction parameters with only a limited subset of interactions requiring quantitative specialization. This framework provides a compact endothelial signaling model for investigating mechanobiological regulation and for future integration with multiscale models of vascular adaptation and disease.

Background: Cerebral autoregulation stabilizes cerebral blood flow over a range of cerebral perfusion pressures, but the precise shape of the pressure-flow relationship remains debated. The classical triphasic pressure-flow relationship was recently challenged by experiments demonstrating a quadriphasic response, hypothesized to arise from vessel-size-dependent pressure-diameter responses. We tested this hypothesis and investigated whether these size-dependent responses originate from heterogeneity in smooth muscle cell (SMC) abundance, SMC behavior, or neither. Methods: We developed a computational multi-scale physics-based model of cerebral autoregulation linking SMC activity to vessel-scale diameter regulation and organ-scale blood flow. Four scenarios were evaluated: passive vessels, homogeneous SMC abundance and behavior, heterogeneous SMC abundance, and heterogeneous SMC behavior. Predicted pressure-diameter responses and pressure-flow relationships were compared across scenarios and against experimental observations. Results: In contrast to passive vessels, homogeneous SMC activation produced partial flow stabilization, highlighting the key role of SMCs in autoregulation. However, only heterogeneous SMC behavior reproduced the experimentally observed vessel-size-dependent trends in pressure-diameter responses. This scenario also showed the best agreement with the experimental organ-scale pressure-flow relationship (R-squared = 0.93, nRMSE = 5.96%). Conclusion: The model suggests that vessel-size-dependent SMC behavior underlies vessel-size-dependent pressure-diameter responses and shapes the relationship between cerebral perfusion pressure and cerebral blood flow.

Reduced-order models are increasingly used to study cardiac physiology and inform patient-specific therapies. However, a model's prediction is only as reliable as its underlying parameters: representative model parameterization is essential to reflect the physiology of the populations these models are meant to represent, including biological sex. Most current models are parameterized from male or sex-agnostic data and/or focus on specific pathologies. Therefore, the goal of this work is to establish a formal parameter estimation pipeline for deriving reduced-order cardiovascular model parameter ranges that are physiologically representative of healthy women and men. We calibrated a closed-loop reduced order model of the heart and circulation separately for healthy female and male populations, using data pooled from eleven healthy cohorts. To account for parameter sensitivity and identifiability, we employed a three-stage parameter subset reduction pipeline: global sensitivity analysis (Sobol's method), collinearity screening (Fisher information matrix), and profile-likelihood identifiability analysis. Sex-specific distributions of parameters that were deemed sensitive and identifiable for each sex, ten for women and 9 for men, were obtained by Hamiltonian Monte Carlo. All the calibrated parameters showed less than 80% overlap between sexes, with the smallest overlap observed in some of the most influential parameters, such as stressed blood volume. Comparing simulations of the calibrated models against allometrically size-matched simulations showed that body size explained some, but not all, of the sex differences. The resulting parameter distributions provide reference ranges usable in future mechanistic and patient-specific simulations to contribute to more inclusive cardiovascular modeling.

The kidney is a critical soft-tissue organ responsible for blood filtration. Accurate constitutive models of the kidney are essential to predict tissue deformation and stress, yet existing models prescribe the strain energy function a priori, largely treat tissue variability as deterministic, and do not distinguish between cortex and medulla. Here we use Gaussian constitutive artificial neural networks to discover probabilistic strain energy functions for kidney cortex and medulla from tension, compression, and shear experiments. The medulla is approximately twice as stiff as the cortex across all three modes: The effective Young's moduli are 3.41, 5.04, and 4.48 kPa for the medulla and 1.43, 2.32, and 2.87 kPa for the cortex in tension, compression, and shear. Both regions display pronounced tension--compression asymmetry. When trained on all three modes, the network discovers strain energy functions in the second invariant alone over the tested deformation range. The cortex model selects two exponential I2 terms, while the medulla model selects an additional linear I2 term. Gaussian external weights propagate specimen-to-specimen variability into closed-form probabilistic stress predictions. Our region-specific probabilistic models capture the mechanical heterogeneity and variability of the kidney and enable more realistic simulations of kidney deformation and stress for surgical planning, needle interventions, and renal trauma.

Tumor spheroid mechanics arise from both the mechanical properties of individual cells and the adhesive interactions that organize them into tissues. The relative contribution of these two factors to the bulk mechanical behavior, however, remains difficult to disentangle experimentally. Here, we develop a computational model of micropipette aspiration to compare the mechanical response of isolated cells and multicellular spheroids within a common computational framework. By independently varying single-cell stiffness and cell-cell adhesion, we quantify their effects on aspiration dynamics, effective elastic modulus, and viscoelastic relaxation. Our results show that increasing single-cell stiffness substantially alters the mechanics of isolated cells but has limited influence on the effective elastic modulus of multicellular spheroids. In contrast, changes in cell-cell adhesion produce pronounced effects on spheroid effective elastic modulus. Nevertheless, both parameters increase the retardation time governing the transition from the initial elastic response to long-time viscous deformation. These findings suggest that multicellular elasticity is governed primarily by intercellular mechanical coupling, whereas the dynamical response to applied stress depends jointly on cell-scale mechanics and cell-cell adhesion.

Cerebral aneurysms are abnormal outpouchings of arteries within the brain and occur in [~]1 in 30 adults. Their initiation, growth, and rupture have been linked to focal blood flow abnormalities--often termed "disturbed" or "hostile" hemodynamics--but commonly-used hemodynamic metrics yield conflicting associations with pathology and lack a unifying mechanistic interpretation. Building on a theoretically-grounded link between wall shear stress and near-wall vorticity, we hypothesized that a topology-based description of near-wall flow can operationalize the concept of hostile hemodynamics in a reproducible way. Inspired by atmospheric tornadic phenomena, we sought a principled taxonomy of coherent near-wall fluid structures with potential mechanobiological and clinical implications. Using high-fidelity computational fluid dynamics simulations in anatomically realistic geometries, we identified coherent near-wall fluid structures whose organization mirrors well-studied atmospheric phenomena: tornado-like columnar rotating cores; downburst-like nonrotating wall-impinging jets with tangential outflow, roll-cloud-like tangential vortices; and mixed configurations. These tornadic events on the aneurysm luminal surface were identified from wall shear stress topology, consistent with its theoretical connection to near-wall vorticity kinematics. The presence of tornadic phenomena--and their imprints on the aneurysm wall--was independently observed in vivo using 4D flow magnetic resonance imaging. By translating concepts from atmospheric physics into vascular biomechanics, this topology-based framework yields a unified mechanistic language for describing near-wall hemodynamics, resolving blood flow complexity into interpretable and reproducible coherent fluid structures, enabling standardized hemodynamic phenotyping, and supporting hypothesis-driven studies of aneurysms and other cardiovascular diseases where greater fluid-mechanical specificity and interpretability may strengthen links between mechanobiology and clinical risk.

Predictive biomechanical models of the tricuspid valve require accurate representation of the chordae tendineae, yet subject-specific chordal geometry is difficult to reconstruct from non-invasive imaging. Here, we adopt a framework for generating functionally equivalent synthetic chordae without prior knowledge of in vivo chordal attachments. To this end, we first adapt an anatomy-informed hyperelastic shape-matching method that establishes correspondence between end-diastolic and end-systolic leaflet configurations using chordal-mimicking forces, a rigid contact template, and a Gaussian-smoothed locally corrective pressure field. We then generate synthetic chordal insertion sites using zone-based rejection sampling and calibrate the unloaded length of each chord by combining reaction forces with the chordal stress-stretch relationship. The framework was evaluated using Texas TriValve 1.1, a high-fidelity finite element model of a human tricuspid valve validated against beating-heart echocardiography. We found that shape matching reproduced the target end-systolic geometry with a mean inter-surface distance of 0.29 {+/-} 0.35 mm. Moreover, synthetic chordal insertions faithfully reproduce end-systolic leaflet deformations. We subsequently examined synthetic chordal configurations containing 202, 225, and 450 insertions, informed by measurements from eight explanted human tricuspid valves. Here, increasing insertion number reduced mean inter-surface distance from 0.63 {+/-} 0.52 mm to 0.49 {+/-} 0.44 mm and reduced contact area errors from 3.52% to 0.59%. Across all configurations, mean maximum principal stretch errors in leaflet belly regions remained below 2.4%, while areal strain errors ranged from 0.34% to 10.01%. These results demonstrate that anatomically informed shape matching coupled with stress-based chordal calibration can reproduce tricuspid valve closure without explicit subject-specific chordal geometry. This framework provides a foundation for generating synthetic subvalvular anatomy for future imaging-derived, predictive tricuspid valve models.

Pulmonary hypertension (PH) is characterized by progressive structural and mechanical remodeling of the pulmonary vasculature, yet few computational frameworks directly link disease mechanisms to longitudinal progression and therapeutic response. In this study, we utilized a multiscale pulmonary arterial growth and remodeling (G&R) framework to capture evolving functional metrics from rat models of PH. This framework couples morphometric tree hemodynamics, constrained mixture theory-based wall mechanics, and maladaptive cellular remodeling. Disease progression was driven by three mechanistically interpretable parameters governing excess smooth muscle production, remodeling activation, and passive stiffening. These parameters were calibrated to longitudinal monocrotaline (MCT) measurements of pressure, wall thickness, and stiffness from prior work using a multiobjective optimization. To show the predictive value of this model, we simulated therapeutic intervention within the same disease-specific framework by using functional cell-level responses to therapy to inform changes in parameter values. Calibration to the study-specific MCT dataset reproduced the temporal increases in pressure, wall thickness, and stiffness, demonstrating that the model could capture multiple features of vascular remodeling simultaneously, with R2 values of 0.81, 0.83, and 0.95, respectively. Simulated treatment reduced pressure, wall thickness, and stiffness. Predicted pressure and wall-thickness responses agreed closely with the corresponding experimental treatment effects, whereas stiffness recovery was overpredicted, suggesting that additional mechanisms may contribute to persistent vascular stiffening after intervention. The framework also captured the overall progression of pulmonary pressure increases across both aggregated MCT and Sugen-hypoxia datasets, suggesting utility across studies and animal models. This work outlines a physics-based, multiscale framework that simulated quantities of direct clinical interest in a mechanistically interpretable platform for linking pulmonary vascular remodeling and treatment response. It supports comparisons across experimental phenotypes and interventions while identifying where constitutive refinements are needed to improve predictive capability across phenotypes.

Cerebrovascular pulsatility drives measurable brain tissue deformation and has been associated with ageing and a range of neurological disorders. Yet how pulsatile haemodynamic forces are transmitted through deformable cerebral arteries into the surrounding brain remains poorly understood, particularly in anatomically realistic vascular geometries. Existing computational approaches have largely treated cerebral fluid and tissue mechanics separately or relied on idealised geometries, limiting our ability to determine how vascular anatomy simultaneously governs intraluminal haemodynamics and extravascular mechanical loading. Here, we develop an image-derived three-dimensional computational framework that jointly resolves pulsatile blood flow, arterial wall deformation and surrounding brain tissue motion in representative cerebral arteries. Four arterial segments, including the middle cerebral artery, middle cerebral artery bifurcation, basilar artery and internal carotid artery, are reconstructed from high-field (5 Tesla) magnetic resonance imaging data of a healthy subject. A finite-deformation fluid-structure interaction model is established by coupling non-Newtonian blood flow, hyperelastic arterial wall and hyper-viscoelastic brain tissue. The predicted tissue response is benchmarked against in vivo magnetic resonance elastography measurements of cardiac-induced volumetric strain over a cardiac cycle. Results reveal spatially localised arterial and tissue deformation whose magnitude and distribution are strongly governed by vascular geometry and wall thickness. Among the segments examined, the internal carotid artery exhibits the largest deformation response, while reduced wall thickness increases strain transmission into the surrounding tissue. Geometrically complex regions also exhibit greater spatial heterogeneity in near-wall haemodynamic metrics. These findings demonstrate that cerebral vascular anatomy simultaneously shapes intraluminal haemodynamics and extravascular mechanical loading. By integrating image-derived vascular anatomy, coupled blood-vessel-brain mechanics and in vivo benchmarking within a unified framework, this study provides a mechanically consistent reference for healthy cerebral pulsatility and establishes a foundation for quantifying how blood-vessel-brain interactions are altered under pathological conditions.

Epithelial monolayers are the foundation of many mammalian organs and have a significant impact on numerous biological processes, including the secretion of cytokines, absorption of waste products and barrier function. During organ development epithelia are subject to many external forces due to the growth of the surrounding tissues. They exhibit complex responses to such forces, with deformation occurring on short timescales (seconds), and relaxation and remodelling occurring on longer timescales (minutes or hours). Existing mathematical and computational models do not typically account for subcellular remodelling, assuming instead that the response to an applied force or strain acts on a single timescale. In this paper we extend an off-lattice, cell-centre modelling framework to account for subcellular and tissue remodelling by introducing a Dynamic Reference Frame (DRF). The DRF provides a discrete, cell-based analogue of morpho-elasticity: it decouples the evolution of each cell's mechanical reference configuration from its instantaneous deformation, providing a phenomenological description of subcellular remodelling processes - cytoskeletal reorganisation, myosin turnover and adhesion bond remodelling - at the cellular scale. Using this extended multicellular model, we reproduce multiscale responses to external mechanical forces, including creep and stress relaxation experiments, and demonstrate that the history of applied deformation influences tissue recovery upon release. Additionally we identify and quantify where and how these multiscale responses occur. The resulting framework allows for more detailed descriptions and analyses of the development and function of biological tissues.

Purpose: Fibrotic remodeling of the uterus, associated with aging, disease, and environmental exposures, alters collagen organization and tissue stiffness, yet how these changes influence organ-level mechanical behavior remains poorly understood. Glutaraldehyde (GA)-induced collagen crosslinking was used as a controlled surrogate for fibrotic remodeling to determine whether image-informed inverse finite element analysis (iFEA), combined with inflation testing and micro-computed tomography (microCT), could detect and quantify the resulting changes in uterine constitutive behavior. Methods: Murine uteri (n = 6 untreated, n = 6 GA-crosslinked) underwent volume-controlled balloon inflation with simultaneous microCT imaging to quantify deformation of the inner and outer wall boundaries for iFEA. Specimen-specific Gasser-Ogden-Holzapfel (GOH) finite element models were optimized by adjusting model parameters to reproduce experimentally measured wall contours throughout inflation. Model performance was evaluated using contour root mean square error (RMSE), and parameter identifiability was assessed through sensitivity analyses. Results: GA treatment significantly increased inflation work, linear stiffness, and maximum inflation resistance (p < 0.001). The iFEA framework accurately reproduced experimental deformation (RMSE < 3%) and revealed significant increases in the estimated GOH parameters C10 (9.2-fold), k1 (2.0-fold), and k2 (2.7-fold), consistent with increased effective tissue stiffness and a shift toward earlier collagen fiber recruitment. Sensitivity analyses demonstrated unique, well-defined minima for all parameter combinations. Conclusion: Image-informed iFEA provides a quantitative framework for relating collagen remodeling to organ-level uterine mechanics through specimen-specific constitutive parameter estimation. This approach establishes a foundation for investigating the mechanical consequences of uterine fibrosis and other remodeling processes.

Background: Ascending thoracic aortic aneurysm (ATAA) is characterized by progressive biomechanical remodeling that is incompletely captured by maximum aortic diameter. The relationships between age, sex, aortic valve phenotype (BAV vs TAV), and aneurysmal mechanics remains incompletely understood. Methods: Ascending aortic specimens from 54 patients undergoing elective ATAA repair were evaluated. Biaxial testing quantified Green strain and local stiffness at membrane tensions of 60 and 120 N/m. Opening angle was used to assess circumferential residual stress and stretch. Associations with age were evaluated using Pearson correlations, and differences between sex and valve phenotype were also explored. Multivariable linear regression models including age, sex, and valve phenotype were used to evaluate their adjusted associations with each mechanical outcome. Results: Increasing age was associated with increased local stiffness and reduced Green strain and remained the strongest determinant of biaxial mechanical properties after adjustment. Compared with tricuspid aortic valve specimens, BAV specimens exhibited lower Green strain together with lower local stiffness. Females demonstrated greater Green strain in unadjusted analyses; however, after adjustment, sex remained associated only with longitudinal Green strain at 60 N/m. None of the evaluated clinical characteristics consistently explained opening angle or residual stretch. Conclusions: Age was the predominant clinical determinant of biaxial mechanical remodeling in ATAA, whereas valve phenotype demonstrated selective associations and sex contributed comparatively little after adjustment. Circumferential residual stress appeared largely independent of these clinical characteristics, suggesting that residual stress and biaxial mechanics capture complementary aspects of aneurysmal remodeling.

Purpose: Atrial fibrosis is associated with thromboembolic risk in atrial fibrillation, but how fibrosis-related contractile impairment alters blood transport in the left atrium (LA) and left atrial appendage (LAA) remains unclear. Prior multi-physics simulations coupled electromechanical (EM) and computational fluid dynamics (CFD) models to show that fibrosis impairs LA function and flow, but were limited to booster-only contraction and could not separate the effects of fibrosis burden from its spatial distribution. Here we examined how fibrosis burden and spatial pattern independently relate to long-term LA and LAA hemodynamics in full-cycle, multi-beat simulations on a common LA anatomy. Methods: Three fibrosis patterns derived from late gadolinium enhancement MRI were considered: one from the patient providing the common anatomy and two from additional patients, with the latter mapped onto the common anatomy via universal atrial coordinates. All three patterns were applied globally at matched burdens of 15.6% and 31.2%, yielding six global-fibrosis cases and a non-fibrotic reference. Two patterns were also confined to the LAA at three local burdens (10%, 25%, and 50%), yielding six LAA-fibrosis cases. A 3D EM model for the LA anatomy was coupled to a 0D circulatory model, so that the simulations resolved the reservoir, conduit, and booster phases. Full-cycle EM-derived wall motion drove subsequent multi-beat CFD simulations. Region-specific emptying fraction (EF), non-dimensional kinetic energy (KE), and the 90th percentile of blood residence time (TR) were evaluated for the LA and LAA. Results: Left ventricular stroke volume remained nearly constant across cases, while increasing global fibrosis burden progressively reduced LA and LAA EF (LA: 36% to 30-31%; LAA: 44% to 28-38%) and monotonically increased LA TR across all patterns. Non-dimensional LA KE showed an approximately quadratic association with LA EF across patterns, whereas LAA KE was more scattered but still increased with LAA EF. In contrast, the relationship between the 90th percentile of LAA TR and LAA EF showed pattern-specific, non-monotonic trends. Conclusions: Fibrosis burden and spatial distribution jointly shape LA and LAA hemodynamics, but through different mechanisms. LAA EF is a robust functional descriptor of chamber-averaged kinetic energy regardless of fibrosis pattern, whereas LAA blood stasis depends additionally on the spatial distribution of fibrosis beyond what bulk emptying fraction predicts. These findings support incorporating imaging-derived fibrosis pattern, and not only burden, alongside functional metrics when assessing thrombogenic risk in the LAA.

Devireddy et al. [1] previously showed that an enthalpy-based, macroscale model simulated tissue freezing histories are in close agreement with a coupled model incorporating cellular water transport and intracellular ice formation (IIF). That comparison, however, was performed for a limited set of tissue geometries, cooling conditions, and biophysical parameters. Here, we revisit that conclusion by reimplementing the published model from its governing equations and systematically expanding the parameter space to determine when microscale processes alter the predicted macroscale freezing response. We incorporate both surface catalyzed nucleation (SCN) and volume catalyzed nucleation (VCN) of intracellular ice with cellular osmotic dehydration and examine a broad range of tissue dimensions, cooling rates, convective boundary conditions, membrane permeabilities, cell sizes, activation energies, intracellular water fractions, and nucleation rates. Linearization of the cellular water transport equation about osmotic equilibrium yields an osmotic relaxation timethat can be compared with the local residence time of tissue undergoing phase change. Their ratio defines a dimensionless group, that organizes the transition between regimes in which the uncoupled enthalpy formulation and the coupled micro and macroscale formulation produce similar or substantially different results. Across the one-dimensional parameter space examined here, the two models remain within a prescribed thermal-error bound for the dimensionless group below approximately 0.1, whereas the discrepancy increases systematically above this threshold. Model divergence also depends on intracellular ice nucleation, i.e., when nucleation is sufficiently slow, cellular water transport can become rate limiting for latent heat release. Conversely, sufficiently rapid intracellular ice formation reduces the macroscopic discrepancy by providing an additional pathway for latent-heat release. Multidimensional simulations further show that geometry alone can move a nominally conventional freezing protocol into the coupled regime. In particular, multidirectional cooling can reduce the local residence time sufficiently to produce substantial differences in predicted frozen volume even when the corresponding thermal histories appear similar. Finally, we formulate an inverse approach for estimating IIF parameters from tissue scale thermal histories. Synthetic simulations indicate that parameter recovery is feasible within a finite range of cooling conditions, providing a testable strategy for future experimental measurements. Together, these results provide a mechanistic criterion for deciding when microscale biophysical processes must be retained in macroscale tissue-freezing models.

Fibrotic strictures are among the most disabling complications of Crohn's disease, permanently narrowing the bowel and impairing motility, yet no approved therapy reverses them. Chronic inflammation alters pacemaker-network coupling, smooth-muscle excitability, and calcium-dependent contractility, while fibrosis thickens the bowel wall, narrows the lumen, and changes tissue mechanics. The relative contributions of these coupled electrical, contractile, and structural alterations to motility loss remain unclear. To address this gap, we develop an integrated electromechanical finite-element framework for fibrostenosing Crohn's disease that couples a fibrosis-driven growth model with a FitzHugh-Nagumo electromechanical model. A full-factorial 25 design of experiments is used to quantify the relative effects of electrical diffusivity, excitation threshold, peak active stress, wall stiffness, and hypertrophic remodeling on cyclic lumen-volume deformation. Motility is quantified by the standard deviation of lumen volume over one contraction cycle. Within the parameter ranges examined, increased wall stiffness emerged as the dominant contributor to motility loss, followed by impaired smooth-muscle contractility. Changes in excitation threshold, hypertrophic remodeling, and electrical diffusivity produced substantially smaller effects. Pairwise interactions were small relative to the dominant main effects, indicating that the mechanisms contributed largely through their individual effects. Our findings suggest that limiting wall stiffening while preserving smooth-muscle contractile function may provide a therapeutic strategy for maintaining intestinal motility in fibrostenosing Crohn's disease.

Large shear deformation in cartilage can lead to tissue failure and injury, which are often precursors to diseases such as post-traumatic osteoarthritis. Under shear, cartilage strain-stiffens prior to failure, which makes stiffening an important indicator of changes to the loading of the extracellular matrix that gives rise to this tissue's properties. While a variety of stiffening mechanisms ranging from uncrimping to fiber-alignment have been explored, whether such mechanisms apply to the depth dependent shear strain stiffening in articular cartilage remains poorly understood. Here, we present experiments and theory that support the idea that the strain-stiffening behavior is caused by a mechanical phase transition from a random isotropic network whose mechanics arise primarily from bending-dominated modes to one where the mechanics arise from stretching-dominated modes. Indeed, consistent with network model predictions, we find that stress-strain data for all depths can be collapsed onto a single universal curve that diverges at a depth-dependent critical strain. These results show that increasing the local concentration of collagen, which constitutes the primary load-bearing network in cartilage, decreases the critical strain at which the strain-stiffening transition occurs, while increasing the low-strain elastic modulus. Furthermore, we demonstrate that bulk compression leads to a higher critical strain and decreased low-strain elastic modulus. Network simulations capture this rich behavior, and indicate that under compression the network collapses, which allows for larger strains before the transition from bending-dominated to stretching-dominated mechanics occurs. Collectively, these findings shed light on the roles that tissue composition and axial strain play in the mechanisms that organize how articular cartilage stiffens during the lead-up to injury.

The mechanical properties of the nucleus are critical for maintaining nuclear integrity and function. We previously showed that chromatin dominates short-extension mechanics whereas lamins provide long-extension strain stiffening. To distinguish the roles of lamin isoforms, micromanipulation nucleus force measurements were performed on isolated nuclei from lamin A/C (Lmna-/-) and lamin B1 (Lmnb1-/-) knockout mouse embryonic fibroblast cells. Lamin A/C knockout does not alter short-extension nuclear stiffness but is essential for strain stiffening at longer extensions. Oppositely, lamin B1 loss reduced short-extension stiffness due to facultative heterochromatin loss while long-extension strain stiffening was slightly increased. Loss of lamin A/C and B1 resulted in similar lamin-chromatin linkers effects as LBR did not change and LAP2{beta} decreased in both. A simulation model of a polymeric lamin shell with stiff lamin A/C and softer lamin B1 subunits can qualitatively recapitulate experimental measurements of lamin knockout cells. Lamin A/C knockout resulted in abnormal nuclear shape but not nuclear blebbing or rupture whereas lamin B1 knockout, similar to other perturbations that cause heterochromatin loss, resulted in increased nuclear blebbing and rupture. This work illuminates the distinct mechanical roles of lamin A/C and B1 in determining nuclear structure and integrity.

Repetitive low-intensity blast (LIB) exposure has been identified as a probable cause of mild blast-induced traumatic brain injury (mbTBI) and a chronic injury risk to U.S. military personnel. However, the human brains biomechanical response to this loading regime remains poorly characterized. Using blast-exposure-validated 3D models of human anatomy, we simulated the intracranial tissue response to blast pressure typically experienced by Warfighters during weapons training. Two scenarios were evaluated, a single-dose exposure and a repetitive-dose exposure, to study intracranial pressure (ICP), shear strains, and spectral content. Ansys LS-DYNA was used to generate planar blast waves with peak overpressures of 4-90 kPa and positive phase durations of 2.2-10 ms. Single exposures produced ICP ranging from 4.7-112.7 kPa, dependent on dose and positive phase duration. Under repetitive LIB exposure, peak ICP increased by 8-26% relative to single exposures, with an increase of high-frequency components (>2 kHz). These results demonstrate that LIB can produce measurable intracranial responses that are amplified through repetition, producing pronounced spectral content and elevated pressures despite low strain levels. This study underscores the need to further investigate cumulative dose effects and the value of computational approaches to clarify hypothesized mbTBI mechanisms in operationally relevant conditions.

Healing of the gastrointestinal (GI) tract following surgical interventions is governed by coupled mechanobiological processes and may lead to severe complications such as perforation or stricture formation. A first comprehensive mechanobiological modeling framework of GI healing is herein proposed. The framework addresses two clinically relevant applications of bioprinting: endoscopic treatment of GI tissue defects and bioprinting-assisted surgical anastomosis. Tissue remodeling is based upon a hyperelastic material model with fiber dispersion coupling the spatiotemporal reaction diffusion dynamics of fibroblasts and cytokines with collagen deposition, reorientation, and plastic deformation. An active stress formulation drives fibroblasts and myofibroblast contraction while a multiplicative decomposition approach is adopted for collagen remodeling as a permanent change. The governing equations are discretised in time with an implicit Euler scheme, finite element in space and solved by a staggered algorithm implemented in the open-source platform FEniCS. Numerical simulations show that the initial collagen fibers content and orientation strongly influence fibroblast infiltration and scar contracture formation. The study quantifies the advantage of having a collagen bioprinted architecture aligned with tissue fibers, to reduce structural remodeling distortion. The model also predicts a strong correlation between collagen deposition and luminal burst pressure during the recovery of mechanical integrity in GI anastomosis.

Ultrasound is commonly used to assess structural changes in symptomatic Achilles tendons, but quantitative biomechanical metrics for progressive tendon deterioration remain limited. The goal of this study was to develop and validate an automated ultrasound tracking algorithm for regional tendon deformation and evaluate strain progression in survived and ruptured tendons during fatigue loading. We hypothesized that maximum strain, average strain, and strain heterogeneity would exhibit different trajectories between groups. Ten cadaveric Achilles tendons underwent cyclic loading with stress tests every 500 cycles until rupture or 150,000 cycles. Ultrasound images acquired during stress tests were analyzed using an automated tracking algorithm to generate spatially resolved regional strain fields. Ultrasound-derived bulk strain was highly correlated with actuator-derived strain in survived (R^2 = 0.968 +/- 0.017) and ruptured tendons (R^2 = 0.972 +/- 0.014). Maximum and average longitudinal strains progressively diverged between groups across fatigue life (Group x FatigueLife: p = 0.003 and p < 0.0001, respectively). During the first 10,000 cycles, average strain decreased in survived tendons ({beta} = -0.0268%, p = 0.0215) but not ruptured tendons ({beta} = 0.0147%, p = 0.1197), with a significant Group x Cycle interaction (p = 0.0061). This study demonstrates that the algorithm quantified Achilles tendon deformation with high fidelity and enabled spatially resolved strain assessment throughout fatigue loading. Maximum and average strain followed different trajectories between groups, whereas strain heterogeneity did not. Early differences in tendon biomechanics suggest that regional strain behavior may change before pronounced differences in absolute magnitude develop.